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Abstract 

We present an exact electrovacuum solution of Einstein-Maxwell equations with infinite sets of 
multipole moments which can be used to describe the exterior gravitational field of a rotating 
charged mass distribution. We show that in the special case of a slowly rotating and slightly 
deformed body, the exterior solution can be matched to an interior solution belonging to the Hartle- 
Thorne family of approximate solutions. To search for exact interior solutions, we propose to use 
the derivatives of the curvature eigenvalues to formulate a C 3 — matching condition from which the 
minimum radius can be derived at which the matching of interior and exterior spacetimes can be 
carried out. We prove the validity of the C 3 — matching in the particular case of a static mass 
with a quadrupole moment. The corresponding interior solution is obtained numerically and the 
matching with the exterior solution gives as a result the minimum radius of the mass configuration. 

Keywords: Matching conditions; Astrophysical Compact Objects. 
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I. INTRODUCTION 



In general terms, the group of massive compact objects is defined in astrophysics to 
contain black holes, neutron stars, white dwarfs, planet-like compact objects, and other 
exotic dense stars. The main common characteristic of compact objects is that they are 
small for their mass. In the case of weak gravitational fields, Newtonian gravity provides an 
adequate physical description of compact objects. However, for the study of very massive 
objects like black holes and neutron stars it is necessary to take into account relativistic 
effects. Moreover, exotic objects like pulsars and quasars can be investigated only in the 
framework of relativistic astrophysics. Furthermore, the accuracy of modern navigation 
techniques made it necessary to take into account relativistic effects of the gravitational 
field of the Earth and other planets of the Solar system. 

In general relativity, the gravitational field of compact objects must be described by a 
metric satisfying Einstein's equations 

iV - \a^R = 8^ , (l) 

in the interior part of the object (T^ u ^ 0) as well as outside in empty space {T^ v = 0). 
Since at the classical level the most important sources of gravity are the mass and charge 
distributions, it is clear that the problem of finding a metric that describes the gravitational 
field of compact objects can be divided into four related problems. The first one consists 
in finding an exact vacuum solution (R^ = 0) that describes the field corresponding to 
the mass distribution. Secondly, the charge distribution must be considered by solving the 
electrovacuum field equations 

R,u = Svr (F^F/ - ±g, u F Xr F x ^ , (2) 

where F^ u is the Faraday tensor corresponding to a charge distribution Q(x^). As for the 
interior region, it is necessary to propose a model to describe the internal gravitational 
structure of the object. The most popular model is that of a perfect fluid 

T^v = (p + p)u^u u - g^p , (3) 

with density p(x^), pressure p(x^), and 4- velocity . Finally, the fourth part of the problem 
consists in considering the mass and charge distribution simultaneously, i.e. one needs an 
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exact solution to the equations 

R^u - ^9^R = 8vr T^(p, p, Q) , (4) 

where a specific model must be proposed for the energy-momentum tensor. Each of the 
above problems is very difficult to solve because Einstein's field equations are highly non- 
linear. To overcome in part these difficulties it is necessary to assume that the gravitational 
field is invariant with respect to certain transformations. In particular, one can assume that 
the field is stationary and axially symmetric, conditions that are in good agreement with 
observations. In fact, compact objects show in general a rotation that is constant over quite 
long periods of time, and their shape deviates only slightly from spherical symmetry. 

In this work, we review a family of electrovacuum solutions that describe the exterior 
gravitational and electromagnetic fields of a rotating compact object. This represents a 
solution to the first two problems mentioned above. It turns out that in the limit of a weak 
gravitational field, a particular solution of this family can be matched smoothly with an 
interior approximate solution. This opens the possibility of searching for an exact interior 
solution that could be matched with the exact exterior solution. As a particular example, 
we consider the case of a static mass with quadrupole moment, and search for a physically 
meaningful interior solution. The problem of matching the solutions arises immediately. It 
is generally accepted that to match an interior and an exterior solution one first should 
determine a matching surface S. If the metrics and their second order derivatives coincide 
on S, the matching is reached in the sense that the curvature of the interior spacetime region 
passes smoothly into the exterior region. This can be considered as a C 2 — matching. In this 
procedure, however, the determination of £ remains unclear. In the investigation of the 
example of an interior solution with quadrupole moment reported in this work, we noticed 
that at certain distance from the origin of coordinates the curvature eigenvalues present a 
non-expected behavior in the sense that they change their sign several times passing through 
points of maximum and minimum local curvature. This behavior could be associated with 
the existence of a repulsive gravitational potential. It then seems plausible that to avoid 
repulsive gravity one can cover the region of repulsion with an interior solution. The location 
of £ can then be anywhere outside the region of repulsion that, in turn, can be invariantly 
defined in terms of the derivatives of the curvature eigenvalues. This can be considered as 
a C 3 — matching procedure. In this work, we use this procedure to match the exterior field 
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of a mass with quadrupole to the corresponding interior solution. 



II. THE EXTERIOR SOLUTION 

Assuming that astrophysical compact objects satisfy the conditions of stationarity and 
axial symmetry, their gravitational and electromagnetic fields can be described by the line 
element |1| 



^2 

ds 2 = f(dt — usdip) 2 — — 



e 2 "(x 2 - y 2 ) ( —j——r + + (* 2 - 1)(1 - y 2 W 



(5) 



x 2 — 1 1 — y 2 

where (t,(p,x,y) are prolate spheroidal coordinates, and the functions /, u and 7 depend 
on x and y only. In general terms, we can say that a solution of the Einstein-Maxwell field 
equations can be considered as a candidate to describe the exterior field of a compact ob- 
ject if the metric functions can be written with the functional dependence / = f(x,y,P n ), 
u = u(x,y,P n ), and 7 = j(x,y,P n ), where P n = (M n , J n ,Q n , H n ), n = 0,1,2,..., repre- 
sent invariant sets of multipole moments associated with all possible sources of gravitation, 
namely, M n are the gravitoelectric, J n the gravitomagnetic, Q n the electric, and H n the 
magnetic moments. A family of solutions satisfying this property was derived in Ref. 2j. 

It was obtained by using the solution generating techniques which are based upon the 
symmetries of the field equations in the Ernst representation 

(£f - TT* ~ 1) V 2 £ = 2(£*V£ - J^VJ 7 ) V£ , (6) 

(Cf - tt* - 1) v 2 ^ = 2(rve - rv^wi 7 (7) 

where the gravitational potential £ and the electromagnetic J 7 Ernst potentials are defined 

as 

1-f-iQ $ 
^ ~ 1 + f + iQ ' ~ 1 + f + iQ ' ( ^ 

with a(x 2 — 1)Q X = f 2 u) y , and a(l — y 2 )Q y = —f 2 u x . Moreover, the potential $ is determined 
uniquely once the electromagnetic potentials A t and A v are given. Finally, V represents the 
gradient operator in prolate spheroidal coordinates. 
The explicit form of the solution can be written as 

(q + + ib + )e 2 ^ + q_ + ib_ 1/2 = ep + zgo 

(a+ + ib + )e 2 ^ - a_ - 1 0) l + £ 
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where 



with 



4> = '£(-l) n q n P n (v)Qn(x) (10) 

n=l 

a± = (x± l) <5_1 [x(l - A/x) ± (1 + A/i)] , (11) 

6 ± = (x±l) 5 - 1 [ 2 /(A + /U ) T (A-/i)] , (12) 

A = ai (x 2 - l) 1 - 5 ^ + J/ ) M - 2 e M E^i(-i)"9nA.- j (13) 

^ = a2 (x 2 - l) 1 ^ - y ) 25 -V*£:° =1 (-i) n *A + , (14) 



and 



^ln^^-gi(x) + P n (y)Q n - 1 (x) 



Pn±= (±l) r 
n-1 

-^(i^^n-fc^tQn-fc+lC^-Qn-fc-l^)] • (15) 

fc=l 

Here P n (y) and Q n (x) represent the Legendre polynomials and functions of second kind, 
respectively. The constant parameters e , go, cr, ati, a 2 , q n , and 5 determine the gravitational 
and electromagnetic multipole moments. The metric functions / and u can be obtained 
from the definitions of the Ernst potentials, whereas the function 7 can be calculated by 
quadratures once / and u are known. In general, this solution is asymptotically flat and 
free of singularities along the axis of symmetry, y = 1, outside certain region situated close 
to the origin of coordinates. The sets of infinite multipole moments can be chosen in such a 
way as to reproduce the shape of ordinary axially symmetric compact objects. 

One of the most interesting solutions contained in this family is the one with non- vanishing 
parameters go = lj Q2 — Q, 5, a\ = a.2 = (cr — m)/a, where m and a are new constants. 
Consequently, the solution possesses the following independent parameters: m, a, 5, and q. 
In the limiting case a = 0, a = 0, q = and 5 = 1, the only independent parameter is m and 
the Ernst potential (Mi determines the Schwarzschild spacetime. Moreover, for a = a = 

h n 

and q = we obtain the Ernst potential of the Zipoy-Voorhees (ZV) [4, |5[ static solution 
which is characterized by the parameters m and 5. Furthermore, for a = a = and 5 = 1, 
the resulting solution coincides with the Erez- Rosen (ER) static spacetime [6]. The Kerr 
metric is also contained as a special case for q = and 5 = 1. The physical significance of 
the parameters entering this particular solution can be established in an invariant manner 
by calculating the relativistic Geroch-Hansen [3, ls| multipole moments. We use here the 



procedure formulated in Ref. js| which allows us to derive the gravitoelectric M n as well 
as the gravitomagnetic J n multipole moments. A lengthly but straightforward calculation 
yields 

M 2k+1 = J 2k = 0, k = 0,1,2,... (16) 
M = m + a(5 - 1) (17) 
M 2 = ^ 3 5q - \cr 3 (5 3 - 35 2 -45 + 6)- ma 2 5{5 - 2) - 3m 2 a(8 - 1) - m 3 , (18) 

J 1 = ma + 2aa(5 - 1) , (19) 
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J 3 = --aa 3 5q 
15 



-a 3 {5 3 - 35 2 - 5 + 3) + ma 2 {35 2 - 65 + 2) + Am 2 a{5 - 1) + m 3 
3 



(20) 



The even gravitomagnetic and the odd gravitoelectric multipoles vanish identically because 
the solution possesses and additional reflection symmetry with respect to the hyperplane 
y = which can be interpreted as the equatorial plane. Higher odd gravitomagnetic and even 
gravitoelectric multipoles can be shown to be linearly dependent since they are completely 
determined in terms of the parameters m, a, q and 5. From the above expressions we 
see that the ZV parameter 5 enters explicitly the value of the total mass M as well as 
the angular momentum Ji of the source. The mass quadrupole M 2 can be interpreted 
as a nonlinear superposition of the quadrupoles corresponding to the ZV, ER and Kerr 
spacetimes. A generalization of the Kerr metric which includes an arbitrary quadrupole 
moment is obtained by imposing the condition 5=1. The resulting multipoles are 

M 2k+1 = J 2k = 0, k = 0,1, 2, ... (21) 

2 / a 2 \ 3/2 

M = m , M 2 = -ma 2 + —qm 3 1 -) , ... (22) 

15 \ m A J 

4 / a 2 \ 3/2 

J\ = ma , J 3 = —ma 3 H qm 3 a 1 , .... (23) 

15 \ m l ) 

It is interesting to note that this particular exact solution in the limit a — > m leads to the 
spacetime of an extreme Kerr black hole, regardless of the value of the quadrupole parameter 
Q- 

In the limiting static case of the ZV metric the only non-vanishing parameters are m = a 
and 5 so that all gravitomagnetic multipoles vanish and we obtain 

M = m5 , M 2 = hm 3 (l - 5 2 ) , (24) 
o 



for the leading gravitoelectric multipoles. Consequently, the ZV solution represents the 
exterior gravitational field of a static deformed body. 

In all the above special solutions the electromagnetic field vanishes identically One can 
easily obtain the corresponding electrovacuum generalizations by assuming that eo 7^ and 
go ^ 0. The computation of the respective electromagnetic multipole moments can be 
performed in an invariant manner and the result can be expressed as 

E n = e M n , H n = g J n . (25) 

This means that the charge distribution resembles the mass distribution. The electric mo- 
ments vanish identically if no mass distribution exists. This result is in accordance with our 
physical intuitive interpretation of a charge distribution. The magnetic moments turn out 
to be proportional to the gravitomagnetic multipoles, with no magnetic monopole. This is 
a physical reasonable result in the sense that the magnetic field is generated by the motion 
of the charge distribution, in the present case, by the rotation of the compact object. 

It is worth noticing that in all the above special solutions we assumed that a\ = a<i 
and obtained generalizations of the Kerr metric with arbitrary quadrupole moment. More 
general solutions can be obtained by relaxing this condition. Consider, for instance, the 
special solution with S — 1, qo — 1, q% — for % — 1, 2, and 

a + ~ a-- V 2 

The resulting potential corresponds to the charged Kerr-Taub-NUT spacetime with total 
charge Q = meo, where / is the Taub-NUT parameter. 

III. MATCHING WITH AN APPROXIMATE INTERIOR SOLUTION 

In general relativity, the gravitational field of a compact object is described by a Rie- 
mannian manifold which must be well defined in the entire spacetime. Usually, the exterior 
metrics are characterized by the presence of curvature singularities in a very specific region 
of spacetime. A possibility to eliminate those singularities is to find a physically meaningful 
regular interior solution that covers the singular region. If the two solutions can me matched 
smoothly on some hypersurface S , one can say that the entire manifold is well defined. In 
this section we will see that a particular solution contained in the above family determines 
a well-defined manifold in the entire spacetime of a particular rotating compact object. 
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Consider the solution (J9]) with the special choice 



e = go = , go = 1 , 9i=0 



q-2 



l-q 



(27) 



i. e., we consider a vacuum solution with a ZV parameter given in terms of the quadrupole 
parameter. Moreover, the angular momentum parameters are chosen as 



cr — m 



«i = a 2 



a" 



2 2 

m — a 



(28) 



Furthermore, let us assume that the condition for a slowly rotating slightly deformed com- 
pact object is satisfied, i. e., 

r,2 



IsjI « kl > 3 = M, 



— ^ « i . 



(29) 



Expanding the resulting potential up to the first order in q and up to the second order in 
a/m, and introducing spherical-like coordinates R = R(x,y) and = Q(x,y), a stationary 
axisymmetric spacetime is obtained whose line element can be written as [10, QJ| 



ds 2 



2M 
~R~ 



1 + 2fciP 2 (cose) + 2 (1 



2.M 



1 J 2 



2A4 \ 



-i r 



1 -2fc 2 P 2 (cose) -2 (1 



R J P 4 
2.M\~ 1 J 2 
~RJ P 4 
J 



(2 cos 2 9-1) J dt 2 
dR 2 



-P 2 [l - 2fc 3 P 2 (cos 6)](de 2 + sin 2 6^ 2 ) + 4^ sin 2 edtdp 

P 



(30) 



where 



J 2 f M\ 5 Q- J 2 /M 
MR 3 V + P / + 8 M 3 



fei- 



6J 2 
P 4 



J 2 5 g - J 2 /M f 
kl + m 4 M 2 R V 



R J \M 



Here Q™ are the associated Legendre functions of the second kind 

Ql(x) = (x 2 - If' 2 
3, o 



3x 2 -2 3 , x + 1 

— 7, x In - 

x 2 — 1 2 x — 1 



Q&x) 



. Sl x + 1 5x — 3x 
1 ) m H — 



2 v y x — 1 x 2 — 1 



(31) 



(32) 



The metric ( 130]) represents the exterior gravitational field of a slowly and rigidly rotating 
compact object in which the deviation from spherical symmetry has been taken into account 
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only up to the first order in the quadrupole parameter q. The main gravitational multipoles 
are the total mass A4, the total angular momentum J and the quadrupole moment Q which 
reduce to 

j2 4 

M. = mil — q), J = —ma, Q = m 3 q . (33) 

m 5 

For the internal structure of the compact object one can assume a perfect fluid model, 
satisfying a one-parameter equation of state, p = p(p). The fact that we are limiting 
ourselves to the case of slow rotation simplifies the problem. In fact, one can first solve the 
non-rotating case and then linearize the field equations with respect to the angular velocity 
Q which is assumed to be uniform. The 4-velocity w M can be identified with the velocity of 
an observer moving with the particles of the fluid so that u R = u @ = 0, and = Qu*. The 
explicit integration of Einstein's equations depends on the particular choice of the density 
function p(R), and on the equation of state. For the sake of simplicity, we consider here 
only the simplest case with p =const so that the pressure p(R) must be determined by 
the field equations. Then, in the limiting case of a non-rotating object the total mass is 
Ai = 47rp7? 3 /3, where 1Z is the radius of the body. Under these conditions, the resulting 
line element can be written as 



ds 2 = (1 + 2$) dt 2 



l + 2i?^P + $ 2 (i?)P 2 (cose) 
dR 



dR 2 



-R 2 [1 + 2$ 2 (i?)P 2 (cose)] [dO 2 + sin 2 Q(d<p - Codt) 2 ], (34) 

where 

$ = $ (i?) + $ 2 (i?)P 2 (cos 6), (35) 

is the interior Newtonian potential. Here $o is the interior Newtonian potential for the 
non-rotating configuration and $ 2 (-R) is the perturbation due to the rotation. Moreover, 
Co = Q — u is the angular velocity of the local inertial frame, where Q =const is the angular 
velocity of the fluid and u = u(R) is the angular velocity of the fluid relative to the inertial 
frame. 

It can be shown that the interior solution (1340 can be matched smoothly on the surface 
R = 1Z with the exterior solution ( l30l) in the special case of an unperturbed configuration 
with $ 2 = and 

MR) = -2vrp (n 2 - y) > ( 36 ) 
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where 1Z is the radius of the non-rotating configuration. In the general case of a rotating 
deformed body, the matching can be performed only numerically. In fact, the solutions 
of the field equations for the inner distribution of mass are calculated using the matching 
conditions as boundary conditions. In this manner, it is possible to find realistic solutions 
that describe the interior and exterior gravitational field of slowly rotating and slightly 
deformed mass distributions. Several examples of this procedure were originally presented 



by Hartle and Thorne 
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13|- 



The above result shows that in the case of a slowly and rigidly rotating compact object 
that slightly deviates from spherical symmetry, the solution presented here is physically 
meaningful and determines a well-defined manifold in the entire spacetime. We interpret 
this result as a strong indication that in general the solution describes the exterior grav- 
itational field of rotating astrophysical compact bodies with arbitrary sets of gravitational 
and electromagnetic multipole moments. 



IV. MATCHING WITH AN EXACT INTERIOR SOLUTION 

Rather few exact stationary solutions that involve a matter distribution in rotation are 
to be found in the literature. In particular, the interior solution for the rotating Kerr 
solution is still unknown. In fact, the quest for a realistic exact solution, representing both 
the interior and exterior gravitational field generated by a self-gravitating axisymmetric 
distribution of a perfect fluid mass in stationary rotation is considered as a major problem 
in general relativity. We believe that the inclusion of a quadrupole in the exterior and in 
the interior solutions adds a new physical degree of freedom that could be used to search 
for realistic interior solutions. We will study in this section the entire Riemannian manifold 
corresponding to the simple case of a static exterior solution with only quadrupole moment. 

The simplest generalization of the Schwarzschild spacetime which includes a quadrupole 
parameter can be obtained from the Zipoy-Voorhees solution with 5 — 1 — q. The corre- 
sponding line element in spherical-like coordinates can be represented as 

ds 2 = (l-^V~V (37) 
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This solution is axially symmetric and reduces to the spherically symmetric Schwarzschild 
metric in the limit q — > 0. It is asymptotically flat for any finite values of the parameters m 
and q. Moreover, in the limiting case m — > it can be shown that the metric is flat. This 
means that, independently of the value of q, there exists a coordinate transformation that 
transforms the resulting metric into the Minkowski solution. From a physical point of view 
this is an important property because it means that the parameter q is related to a genuine 
mass distribution, i.e., there is no quadrupole moment without mass. To see this explicitly, 
we calculate the multipole moments of the solution by using the invariant definition proposed 
by Geroch [7|. The lowest mass multipole moments M n , n — 0, 1, . . . are given by 



whereas higher moments are proportional to mq and can be completely rewritten in terms 
of Mo and This means that the arbitrary parameters m and q determine the mass 
and quadrupole which are the only independent multipole moments of the solution. In the 
limiting case q = only the monopole M = m survives, as in the Schwarzschild spacetime. 
In the limit m = 0, with q 7^ 0, all moments vanish identically, implying that no mass 
distribution is present and the spacetime must be flat. This is in accordance with the result 
mentioned above for the metric f !38p . Furthermore, notice that all odd multipole moments 
are zero because the solution possesses an additional reflection symmetry with respect to 
the equatorial plane. 

We conclude that the above metric describes the exterior gravitational field of a static 
deformed mass. The deformation is described by the quadrupole moment M 2 which is 
positive for a prolate mass distribution and negative for an oblate one. Notice that in order 
to avoid the appearance of a negative total mass M the condition q < 1 must be satisfied . 

A. Matching conditions 

In this subsection we analyze several approaches which could be used to determine the 
matching hypersurface S. Instead of presenting a rigorous analysis, we will present an 
intuitive method based on the behavior of the curvature and the motion of test particles. 

To investigate the structure of possible curvature singularities, we consider the 




(38) 
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Kretschmann scalar K = R^ v \ t R^ vXt . A straightforward computation leads to 



K 



16m 2 (1 - q) 2 (r 2 - 2mr + m 2 sin 2 8) 2q2 - 4q - 
r 4(2-2 9 + 9 2) (1 - 2m/r) 2 (« 2 -« +1 ) 



L(r,0) , 



(39) 



with 



L(r,8) 



3(r — 2m + qm) 2 {r 2 — 2mr + m 2 sin 2 8) 

—q{2 — g) sin 2 8[q 2 — 2q + 3(r — m)(r — 2m + gm)] . 



(40) 



In the limiting case q = 0, we obtain the Schwarzschild value K = 48m 2 /r 6 with the only 
singularity situated at the origin of coordinates r — > 0. In general, one can show that 
the singularity at the origin, r = 0, is present for any values of q. Moreover, an additional 
singularity appears at the radius r = 2m which, according to the metric ( )38|) . is also a horizon 
in the sense that the norm of the timelike Killing tensor vanishes at that radius. Outside the 
hypersurface r = 2m no additional horizon exists, indicating that the singularities situated 
at the origin and at r = 2m are naked. Moreover, for values of the quadrupole parameter 
within the interval 



from the origin of coordinates. This type of singularity is always contained within the naked 
singularity situated at the radius r = 2m, and is related to a negative total mass Mq for 
q > 1. Nevertheless, in the interval q e (1 — ^3/2, 1]\{0} the singularity is generated 
by a more realistic source with positive mass. This configuration of naked singularities is 
schematically illustrated in Fig. [TJ 

The analysis of singularities is important to determine the matching hypersurface S. 
Indeed, in the case under consideration it is clear that S cannot be situated inside the 
sphere defined by the radius r = 2m. To eliminate all the singularities it is necessary to 
match the above solution ( 138]) with an interior solution which covers completely the naked 
hypersurface r = 2m. 

Another important aspect related to the presence of naked singularities is the problem 
of repulsive gravity. In fact, it now seems to be established that naked singularities can 




(41) 



a singular hypersurface appears at a distance 




(42) 
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FIG. 1: Structure of naked singularities of a spacetime with quadrupole parameter q. Plot (a) 
represents the limiting case of a Schwarzschild spacetime (q = 0) with a singularity at the origin of 
coordinates surrounded by the horizon (dashed curve) situated at r = 2m. Once the quadrupole 
parameter q is included, the horizon transforms into a naked singularity (solid curve) and the 
central singularity becomes naked as well. This case is illustrated in plot (b). For values of the 
quadrupole parameter within the interval q £ (l - v/3/2, 1 + v/3/2) \{0}, two additional naked 
singularities appear as depicted in plot (c). 

appear as the result of a realistic gravitational collapse [ijj and that naked singularities can 
generate repulsive gravity. Currently, there is no invariant definition of repulsive gravity in 



the context of general relativity, although some attempts 



lave been made by using invariant 



15Ml71]. Nevertheless, it is possible to 



quantities constructed with the curvature of spacetime [1 
consider an intuitive approach by using the fact that the motion of test particles in stationary 
axisymmetric gravitational fields reduces to the motion in an effective potential. This is a 
consequence of the fact that the geodesic equations possess two first integrals associated 
with stationarity and axial symmetry. The explicit form of the effective potential depends 
also on the type of motion under consideration. 

In the case of a massive test particle moving along a geodesic contained in the equatorial 
plane (9 = n/2) of the Zipoy-Voorhees spacetime (138|) . one can show that the effective 
potential reduces to 



V eff ~ [ 1 r 



i4('--r 

r z V r J 



(43) 



where L is constant associated to the angular momentum of the test particle as measured by 
a static observer at rest at infinity. This expression shows that the behavior of the effective 
potential strongly depends on the value of the quadrupole parameter q. This behavior is 
illustrated in Fig. |2j 
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FIG. 2: The effective potential for the motion of timelike particles. Plot (a) shows the typical 
behavior of the effective potential of a black hole configuration with q = 0. The case of a naked 
singularity with q = 1/2 is depicted in plot (b). 

Whereas the effective potential of a black corresponds to the typical potential of an at- 
tractive field, the effective potential of a naked singularity is characterized by the presence 
of a barrier which acts on test particles as a source of repulsive gravity. Although this re- 
sult is very intuitive, the disadvantage of this analysis is that it is not invariant. In fact, a 
coordinate transformation can be used to arbitrarily change the position of the barrier of 
repulsive gravity. Moreover, the identification of the spatial coordinate r as a radial coor- 
dinate presents certain problems in the case of metrics with quadrupole moments jl8 |. To 
avoid this problem we investigate a set of scalars that can be constructed from the curvature 
tensor and are linear in the parameters that enter the metric, namely, the eigenvalues of the 
Riemann tensor. Let us recall that the curvature of the Zipoy-Voorhees metric belongs to 
type I in Petrov's classification. On the other hand, type I metrics possess three different 
curvature eigenvalues whose real parts are scalars [49j] . The explicit calculation of the curva- 
ture eigenvalues for this metric shows [20] that all of them are real and, consequently, they 
behave as scalars under arbitrary diffeomorphisms. The resulting analytic expressions are 
rather cumbersome. For this reason we performed a numerical analysis and found out the 
main differences between black holes and naked singularities. The results are illustrated in 
Fig. El 

We took a particular eigenvalue which represents the qualitative behavior of all the eigen- 
values. In the case of a black hole, the eigenvalue diverges near the origin of coordinates, 
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FIG. 3: Behavior of the curvature eigenvalue on the equatorial plane (9 = vr/2) of the Zipoy- 
Voorhees metric. Plot (a) corresponds to a black hole solution with q = 0. Plot (b) illustrates the 
behavior in case of a naked singularity with q = —2. 

where the curvature singularity is situated, and it decreases rapidly as r increases, tending 
to zero at spatial infinity. In the case of a naked singularity the situation changes drastically. 
The eigenvalue vanishes at spatial infinity and then increases as the value of the radial coor- 
dinate decreases. At a specific radius r = r min , the eigenvalue reaches a local maximum and 
then rapidly decreases until it vanishes. This oscillatory behavior becomes more frequent as 
the origin of coordinates is approached. It seems plausible to interpret this peculiar behavior 
as an invariant manifestation of the presence of repulsive gravity. On the other hand, if one 
would like to avoid the effects of repulsive gravity, one would propose r m j ra as the minimum 
radius where the matching with an interior solution should be carried out. If we denote the 
eigenvalue as A, then r min can be denned invariantly by means of the equation 



Or 



. (44) 



Then, the radius r m j„ determines the matching hypersurface E and one could interpret 
condition (|44l) as a C 3 — matching condition. In concrete cases, one must calculate all possible 
eigenvalues Aj and all possible points satisfying the matching condition dXi/dr = 0. The 
radius r m i n corresponds then to the first extremum that can be found when approaching the 
origin of coordinates from infinity. In the next section we will show that this approach can 
be successfully carried out in the case of the Zipoy-Voorhees metric. 
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B. An interior solution 



In the search for an interior solution that could be matched to the exterior solution with 
quadrupole moment given in Eq. (|38|) . we found that an appropriate form of the line element 
can be written as 

P 2 7 ( a t i \ ,,2 

ds 2 = fdt 2 - e — I ^- + d6 2 J - fjdtf? , (45) 

where 

e 2 7 o = ( r 2 _ 2mr + m 2 cog 2 d yy(r,9) ? ^ 

and / = f(r, 6), h = h(r), and fi = fi(r, 9). This line element preserves axial symmetry and 
staticity. 

The inner structure of the mass distribution with a quadrupole moment can be described 
by a perfect fluid energy-momentum tensor ([3]). In general, in order to solve Einstein's 
equations completely, pressure and energy must be functions of the coordinates r and 8. 
However, if we assume that p = const, the resulting system of differential equations is still 
compatible. The assumption of constant density drastically reduces the complexity of the 
problem. Then, the corresponding field equations reduce to 

Pr = -^(P + P)j, P9 = ~\(P + P)j > ( 47 ) 

1 / /ze 270 \ 

^ rr = ~2h \^ ee + hr ^ r ~ 327rp— J~ J ' ( 48 ) 

f ft IK Hr\ ft pefe fee . Q „ (3p + p)e 2 ^ 

frr = 7-{2h + j) fr+ hf-^h-T + 8n h • (49) 

Moreover, the function 7 turns out to be determined by a set of two partial differential 
equations which can be integrated by quadratures once / and \x are known. The integrability 
condition of these partial differential equations turns out to be satisfied identically by virtue 
of the remaining field equations. 

Although we have imposed several physical conditions which simplify the form of the field 
equations, we were unable to find analytic solutions. However, it is possible to perform a 
numerical integration by imposing appropriate initial conditions. In particular, we demand 
that the metric functions and the pressure are finite at the axis. Then, it is possible to plot 
all the metric functions and thermodynamic variables. In particular, the pressure behaves 
as shown in FigJU 
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FIG. 4: Plot of the inner pressure as a function of the spatial coordinates. 

It can be seen that the pressure is finite in the entire interior domain, and tends to zero 
at certain hypersurface R(r, 9) which depends on the initial value of the pressure on the axis. 
Incidentally, it turns out that by increasing the value of the pressure on the axis, the "radius 
fuction" R(r, 8) can be reduced. Furthermore, if we demand that the hypersurface R(r, 6) 
coincides with the origin of coordinates, the value of the pressure at that point diverges. 
From a physical point of view, this is exactly the behavior that is expected from a physically 
meaningful pressure function. 

This solution can be used to calculate numerically the corresponding Riemann tensor and 
its eigenvalues. As a result we obtain that the solution is free of singularities in the entire 
region contained within the radius function R(r,9). In particular, one of the eigenvalues 
presents on the equatorial plane the behavior depicted in Figj5j All the eigenvalues have a 
finite value at the symmetry axis and decrease as the boundary surface is approached. 

To apply the C 3 — matching procedure proposed above we compare the behavior of the 
eigenvalue plotted in FigfS] with the corresponding eigenvalue plotted in Figj5j using the 
same scale in both graphics. The result is illustrated in FigJHJ It then becomes clear that 
the first possible point where the matching can be performed is exactly at r m i n which in 
this particular case corresponds to r m j n w 5Mo. This fixes the initial value of the pressure 
on the axis which is then used to attack the problem of matching the interior and exterior 
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FIG. 5: Behavior of the curvature eigenvalue on the equatorial plane (6 = tt/2) of the interior 
solution. 




FIG. 6: Curvature eigenvalues of the interior solution and of the exterior solution with the same 
scale. 

metric functions. In all the cases we analyzed, we obtained a reasonable matching, withing 
the accuracy of the numerical calculations. We repeated the same procedure for different 
values of the angular coordinate (9 = n/A and = 0), and obtained that the matching can 
always be reached by fixing in an appropriate manner the arbitrary constants that enter the 
metric functions / and ji. 
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V. CONCLUSIONS 



In this work, we presented an exact electrovacuum solution of Einstein-Maxwell equations 
which contains four different sets of multipole moments. An invariant calculation shows 
that they can be interpreted as the gravitoelectric, gravitomagnetic, electric and magnetic 
multipole moments. The solution is asymptotically flat and is free of singularities in a 
region situated around the origin of coordinates. The rotating Kerr metric is contained as a 
special case. The NUT parameter can also be included by a suitable choice of the arbitrary 
constants which enter the Ernst potentials. We conclude that this solution can be used to 
describe the exterior gravitational field a charged rotating mass distribution. 

In the particular case of slowly rotating and slightly deformed mass distribution we ob- 
tained the explicit form of the metric, and showed that it can be matched with an interior 
solution which is contained within the class of Hartle-Thorne solutions. This reinforces the 
conclusion that the solution represents the interior as well as the exterior gravitational field 
of astrophysical compact objects. 

We study the problem of matching the interior and exterior spacetimes. We propose a 
C 3 — matching which consists in demanding that the derivatives of a particular curvature 
eigenvalue are smooth on the matching hypersurface. To prove the validity of this approach 
we derived an interior solution for the simplest case of a static mass with an arbitrary 
quadrupole moment, represented by the Zipoy-Voorhees vacuum solution. The numerical 
integration of the corresponding field equations shows that interior perfect fluid solutions 
exist which are characterized by a constant density profile with a variable pressure. Fixing 
the value of the angular coordinate 9, we performed numerically the C 3 — matching. As a 
result we obtain a minimum radius at which the matching can be carried out and a fixed 
value for the pressure on the symmetry axis. These values are then used to reach the 
smooth matching of the interior and exterior metric functions. In all the cases analyzed in 
this manner we obtained a reasonable numerical matching. 

The idea of using the C 3 — matching condition to determine the minimum radius, at 
which an interior solution can be matched with an exterior one, has been proved also in a 
particular case where analytical methods can be applied, namely, in the case of the Kerr- 
Newman class of solutions. The obtained results are reasonable and compatible with other 
results obtained by analyzing the motion of test particles 21|. These results indicate that 
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it should be possible to determine the minimum radius of an astrophysical compact object 
by using the idea of the C 3 — matching presented here. To prove this conjecture in general, 
it will be necessary to use more powerful methods related to the mathematical behavior of 
geodesies and curvature. This problem is currently under investigation [22|. An important 
application of this analysis would be to relate the minimum size of a compact object with 
its binding energy. As a result we would obtain the maximum binding energy which is 
physically allowed for an astrophysical compact object. 
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